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Abstract 

Ab initio approaches are introduced for calculations of inelastic losses and vibrational damping in 
core level x-ray and electron spectroscopies. From the dielectric response function we obtain system- 
dependent self-energies, inelastic mean free paths, and losses due to multiple-electron excitations, 
while from the dynamical matrix we obtain phonon spectra and Debye- Waller factors. These 
developments yield various spectra and optical constants from the UV to x-ray energies in aperiodic 
materials, and significantly improve both the near edge and extended fine structure. 
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Theories of x-ray and e 



creasingly sophisticated 



^ctron spectra and other "optical constants" have become in- 
For example, treatments of such spectra can now go beyond 
the independent-particle approximation, e.g., with time-dependent density functional theory 
(TDDFT) or the Bethe-Salpeter equation (BSE) P, 0]. However, less attention has been 
devoted to the effects of inelastic losses and vibrational damping, especially at high energies. 
These losses are inrportant b ot h in Sp ectro SCOpy and other a Ppll cat 10ns QQ. T huS current 
treatments 

Baa 

often utilize simplified models, e.g., semi-empirical or electron-gas self- 
energies and Einstein or Debye models for vibrations, with considerable variations in the 
results. 

To address these difficulties we introduce in this Letter ab initio approaches for these 
many-body damping effects. We first present an approach for calculating the dielectric re- 
sponse and optical properties within a real-space Green's function (RSGF) formalism 0,0]. 
From these results, we derive a many-pole model for the dielectric function, which is then 
used to obtain photoelectron self-energies, inelastic mean free paths (IMFPs), and contribu- 
tions from multi-electron excitations. Next, from calculations of the dynamical matrix we 
obtain phonon spectra and quantitative Debye- Waller (DW) factors. These developments 
enable a number of applications, e.g., improved calculations of spectra and inelastic losses 
from the UV to hard-x-ray limits. 

A key physical quantity in these calculations is the dielectric response function x( u ) Cl- 
in the long wave-length limit, various spectra are related to x{ u ) through the local atomic 
polarizability a(u) [9], 

a(u>) = / dxdx' d){r) %(r, r', u) d(r) 

(1) 

= Mi,L (G l ,l'(E)) Mjji . 

i,L,L' 

Here, M^l are screened dipole matrix elements between occupied core states \i) and final 
scattering states \L, 0), and Gl,l'{E) are matrix elements of the photoelectron Green's func- 
tion in an ang * nronrentntn' and site oa Sls Q. Tire o P erator d W re« tire drpoie 
coupling to the photon field, and the brackets (• • • ) denote a configurational average, which 
can be expressed in terms of DW factors. From a(u) we obtain the dielectric constant 



e{uj) = 1 + 47ma(u), where n = N/V is the atomic number density 12]. Other optical con 



stants (e.g., absorption coefficient /i, reflectivity 71, anomalous x-ray scattering amplitudes 
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FIG. 1: Energy loss function — Imje ] for foe Cu vs photon energy u from this work (solid), and 



from experiment (a) (dash-dot) 



10| and (b) (dash-dot-dot) 



fi + if2, etc.) are related to e = e\ + ie-i- Typical results, e.g., for the energy loss function 
— Im [e _1 ] = e 2 /[e^ + e|] for Cu (Fig. UJ), are in good agreement with experiment 10, [h| in 
the UV and beyond. This and other examples based on our approach are tabulated on the 

www 0- 

The quantities in Eq. ((TJ) are calculated using an extension of our RSGF code [3] that 
includes corrections to the independent-particle approximation [4j. Compared to other ap- 
proaches, the RSGF method is advantageous for core-level spectra since it avoids an ex- 
plicit calculation of energy eigenstates and is applicable to periodic and aperiodic systems 
alike. For each core state \i), we first calculate the atomic-like background with no scat- 
tering; fine structure is then added with full multiple scattering (MS) to about 50 eV, and 
then with the MS path expansion to about 1500 eV above each edge. In our code both 
G(E) = 1/(E — hi + i0 + ) and the scattering states \L, 0) = Ri{t)Yl{t) are obtained with a 
self-consistent final-state, muffin-tin Hamiltonian h! = p 2 /2 + V coul + H{E). This Hamilto- 
nian includes a screened core-hole, an energy-dependent self energy and the core- hole 

lifetime T (we use Hartree atomic units e = H = m = 1). Our screened core-hole agrees 
well with that from the random phase approximation, and gives results for core-level spectra 
which are a good approximation to the BSE 14]. The screened dipole operators in M ac- 
count for the induced local dipole fields jj. For metals, the low-frequency behavior of e{uS) 
is approximated by an additive Drude term; otherwise there are no adjustable parameters. 

From the dielectric response, we obtain the self-energy E(£") using the "GW approxima- 
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tion" |1_ 

E(E) = i J ^-G(E - u)W(u;)e-^, (2) 

where W = e~ 1 (c<j)l / is the screened coulomb interaction, and matrix indices (r, r') are 
suppressed. For G(E), we use the free Green's function, ignoring MS terms. Hence our 
represents an average self-energy, which is nevertheless adequate for many applica- 
tions. Current GW implementations for core spectra P] typically use a single-plasmon pole 
approximation for the dielectric function. Although computationally efficient, Fig. 1 shows 
that a single pole may not be a good representation. Thus, to obtain improved self-energies 
that exploit the efficiency of pole models, we match our ab initio e(uj) to a many-pole model, 
with poles spread over a broad spectral range, 

-Im [e -1 (q,o;)] = ^J2 9 ^ 5 ( U ~ u 3 ( q ))" ( 3 ) 

3 

Here gj = (2Aujj/7iujj) Im [e _1 (uJj)] is the strength of pole j and Auj is the pole spacing. 
For simplicity, the momentum dependence of each pole is approximated by a polynomial in 
q 2 : u}j(q) = \ui 2 + q 2 /3 + q 4 /4l ^ 2 [l6| . This approximation is roughly consistent with explicit 
calculations 1171 . but the precise dispersion is not critical since q and uj are integrated 

n n n 

over [16[. Many-pole models [18] and other methods [19( have been developed for more 
accurate self-energy calculations; however, these methods are computationally demanding 
and not applicable at high photoelect™ eae rgl e s (e.g., fa the extended x- r ay absorption fine 
structure (EXAFS)). Our model is similar to that of Ref. |7J, but does not rely on empirical 
optical constants. Thus using Eq. (J3J) in Eq. (J2J), we obtain an efficient approximation for 
the self-energy as a sum over single-pole models with excitation energies Uj, i.e., = 
Y^j9j ^(E.Uj). As an application we have calculated the IMFP A = (E /2) 1 / 2 /|Im£(£)| 
for various materials Q|. Fig. El shows A for fee Cu using a 100-pole model matched to the 
loss function in Fig. ^ Clearly this model corrects the excessive loss of the single-pole model 
below 100 eV, and yields quantitative agreement with semi-empirical fits to experiment 
5], |^ . Our self-energy is also in reasonable agreement with that of Ref. . A similar 



approach 



2l| can be applied to inelastic electron scattering (f| and to photoemission spectra 



Our many-pole model also permits calculations of inelastic losses due to multi-electron 
(e.g. "shake-up" and "shake-off") excitations. These excitations correspond to satellites 
beyond the quasi-particle peak in the spectral function A = (— ImG cff . Here G e s 
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FIG. 2: Inelastic mean free path A for fee Cu vs electron wave number k = (2-E 1 ) 1 / 2 and kinetic 
energy E for our many-pole model (solid), a single plasmon-pole model (dashes); and fits to 
experiment (dash-dots) (a) Q and (b) Q. 

is an effective, one-particle propagator, which can be calculated using a generalization of 
the GW approximation with similar ingredients ,22]. In addition to intrinsic losses, G c g 
contains energ y-d ependent interference terms which tend to suppress the satellites. As 
shown in Ref. |22j], these excitations can be included in /i(cu) in terms of a convolution of 
a normalized spectral function A and the quasi-particle absorption coefficient fi qp in the 
absence of satellites, 



du'A(u,uj')fi qp (uj - to') = (fi qp (uj)), 



(4) 



where u' is the excitation energy. Likewise, the net EXAFS is given by a convolution of A 
with the quasi-particle fine structure. Thus for each path R in the MS path expansion, the 
convolution over the oscillatory fine structure yields an amplitude reduction factor Sj^u) 
and (since A is asymmetric) a negative phase shift Sr(u), 



<e 2 ^> = S» e 



2ikR+5 R (uj) 



(5) 



For the first shell of Cu, we find that S R (u) and Sr(lv) cross-over smoothly from the adi- 
abatic (or quasi-particle) limit at threshold to nearly constant values (~0.92 and ~ —0.2 
rad, respectively) in the EXAFS, consistent with experiment S R ~ 0.9 (3]. Clearly, our 
treatment of losses yields significant improvements in amplitude and phase compared to the 
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FIG. 3: X-ray absorption fx vs photon energy to for fee Cu for our many-pole model with inelastic 
losses (solid); a single-pole model (dashes); and experiment (dash-dots) |23j |. 
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FIG. 4: Fourier transform of the Cu K-edge EXAFS x(i?) at 10K vs half MS path-length R with 
(solid), and without inelastic losses (dashes), and from experiment (dash-dots) 2^ |. 

single-pole model, both for the near edge spectra (XANES) (Fig.EJ) and the EXAFS (Fig.HJ). 
Additional details are given elsewhere 21]. 

Finally, we present an ab initio approach for calculating the effects of vibrational damping 
on the spectra. Our method is based on calculations of the dynamical matrix D together with 
a Lanczos algorithm for the phonon spectra (3, Q • Instead of relying on spring-constants 
or semi-empirical approximations, we use first principles electronic structure methods to 
determine the real space matrix elements of D, 

d 2 U 



dUjladUj'l 



■j'l'p 



(6) 



where U is the total internal energy of the system and Uji a is the a = {x, y, z} Cartesian 
displacement of atom j of mass Mj in unit cell /. The DW factor exp (—2a R k 2 ) for a given 
MS path R is related to the the second cumulant a 2 R = ((5R) 2 ) of the vibrational distribution 
function for that path 24^ . 

°1{T) = - pn^coth^du. (7) 
Pr Jo * 

Here p R is a reduced mass for path R, (3 — l/k B T and p R {u 2 ) = (Q R \S {oj 2 1 — D)| Qr) 
is the projected density of phonon modes, which is approximated by a six-step Lanczos 
recursion, starting from a path displacement seed state \Qr)- Other phonon spectra can be 
obtained by varying the seed Although our approach is applicable both to molecules 
and condensed systems, we present results here only for two crystalline examples. For 
these systems the Dji a yiip are obtained using ABINIT [27j . Briefly, the reciprocal-space 
dynamical matrix Djaj'p (q) j3| is calculated on a 4 x 4 x 4 grid of q-points and smoothly 
interpolated within the Brillouin zone. The real-space force constants (i.e., second derivatives 
of U) are then obtained from an inverse Fourier transform. These calculations use Troullier- 
Martins-FHI pseudopotentials 29] with energy cutoffs of 60 and 15 Hartrees for Cu and Ge 
respectively, while an 8 x 8 x 8 Monkhorst-Pack grid is used for the electronic energies. Table 
presents results obtained with local density approximation (LDA) 3l| and generalized 
gradient approximation (GGA) 3^ functionals. The results all lie within about ±10% 
of experiment 3^|, with the LDA usually closer. For comparison we show results for the 
correlated Debye (CD) model 2^, which works well for Cu, but gives significant errors for 
anisotropic structures like Ge. We also calculate anharmonic corrections. The first cumulant 
is obtained from the net thermal expansion yiel ding 0.0217 A for fee Cu at 300K, in accord 
with the experimental value 0.0205 ± 0.0009 A 3J|. The third cumulant is then estimated 
from relations among the cumulants j^]. Further details are given elsewhere 

In summary, we have developed efficient, ab initio approaches for calculations of the key 
many-body damping factors in core-level x-ray and electron spectroscopies. Our calculations 
of inelastic losses are based on a many-pole representation of the dielectric function and an 
extension of our RSGF code that includes corrections to the independent-particle approx- 
imation. This approach yields significantly improved self-energies compared to single-pole 
models, as well as quantitative IMFPs and estimates of losses due to multi-electron exci- 
tations. We have also developed an efficient approach for calculations of phonon spectra 
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TABLE I: Debye- Waller factors (in 1CT 3 A 2 ) for fee Cu and dia Ge at 295K for leading near- 
neighbor shells. 



Cu 



Shell 


LDA 


GGA 


CD 


Exp. [30] 


1 


8.55 


7.15 


8.91 


7.93 ±0.16 


2 


11.06 


9.50 


10.99 


11.08 ±0.50 


3 


10.80 


9.14 


11.32 


9.58 ±0.60 






Ge 






Shell 


LDA 


GGA 


CD 


Exp. [30] 


1 


3.77 


3.41 


5.12 


3.5 ±0.15 


2 


10.21 


11.59 


7.43 


9.3 ±1.1 



and MS DW factors. All of these many-body effects are important in applications to op- 
tical constants over a broad spectrum, and yield improved amplitudes and phases from 
the UV to x-ray energies [13]. Since our approach includes solid state effects (e.g., edge 
shifts, fine structure, and temperature dependent D W factors) and is applicable to general 
aperiodic materials, these results complement and can potentially replace empirical tables 



aperiodic materials 

h h y , y, q 



or atomic models 



38| for many applications. Extensions, e.g, to spec- 



tra at finite momentum transfer and to the visible regime with improved potentials are in 
progress. 
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